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ABSTRACT 

We set up equations and boundary conditions for magnetic rotator winds and 
^ _ disks in axially symmetric hot stars in a steady state. We estabhsh a theorem 

' stating that if a disk region has no meridional motion but its magnetic field has 

a normal component at a point Q on its shock boundary, the angular velocity 
>• ■ of the disk region at Q is the same as the angular velocity of the star at the 

point Pj, at which the magnetic field line through Q is anchored. When there 
' is equatorial symmetry, all points of the disk along the field line through Q will 

Q ■ have the same angular velocity as P^. Also, we show that for a given value of 
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the magnetic field strength, if the rotation rate is too high or the flow velocity 
into the shock boundary is too low, a Keplerian disk region will not be formed. 



Qh' We consider the formation of disks in magnetic rotators through the processes 

O ' of fill-up and diffusion into Keplerian orbits. At the end of the fill-up stage 

• the density of the disk increases significantly and the magnetic force in the disk 

■ becomes negligible. If the meridional component Bm of the field at the surface 

is larger than a minimum value Bm,min, a Keplerian disk can form. The radial 
^ ' extent of the Keplerian region is larger for larger values of Bm and is largest 

o3 ■ when Bm equals an optimal value Bm,opt- The extent does not increase when Bm 

is larger than Bm,opt- If a is the ratio of rotational speed to the critical rotation 
speed at the photosphere, the inner and outer radii of the maximal quasi-steady 
Keplerian disk region are given by a~^/^P and 2^^^a~'^^^R, respectively, where R 
is the stellar radius. For models with dipole-type fields, the values of Bm,min in 
B-type stars are of order 1 G to 10 G and in 0-type stars they are about 500 G. 
Because the values of Bm required for disk formation in B-type stars are relatively 
small and the fill-up time is short, we suggest that meridional circulation may 
play a role in some of the time-variation observed in disks of Be stars through 
its effect on the magnetic field near the photosphere. 



Subject headings: stars: winds, outflows — stars: magnetic fields — stars: rota- 
tion — stars: emission-line. Be — circumstellar matter 
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1. Introduction 

This paper has three main objectives. The first is to develop a framework of equations 
and boundary conditions with general validity for magnetic rotator winds and disks in hot 
stars. The second is to determine the conditions under which Keplerian disk regions may be 
formed. The third is to develop models to explain the formation and properties of Keplerian 
disk regions. In a future paper we intend to consider applications to Wind Compressed 
Disk(WCD) regions. Earher studies of magnetic rotator winds, such as the one due to Weber 
& Davis (1967), focussed on one-dimensional models involving equatorial winds mainly to 
consider magnetic braking. Mestel (1968) set up a system of equations for axially symmetric 
stars with dipole fields. His equations were based on the solutions obtained by Chandrasekhar 
(1956) for steady state magnetohydrodynamic (MHD) systems. Sakurai (1985) developed 
a numerical method for two-dimensional magnetic rotator winds using the Weber-Davis 
approach. There have been several other studies of equatorial and non-equatorial winds and 
flows in the context of disk formation, and a review of such work is given in a recent paper by 
Cassinelli et al (2002) in which a Magnetically Torqued Disk (MTD) model was proposed to 
explain the Be star phenomenon. An important difference between the MTD model and our 
models is that we consider weaker magnetic fields that are sufficient to provide the torque 
required to sustain higher rotation rates in the fiows towards the disk region, whereas the 
fields discussed in the MTD paper are strong fields required to spin up the disk. Since the 
disk material is much denser than the the pre-shock fiow material, the fields required in the 
MTD model are stronger than the fields required in our models. 

Over the past decade, several models have been developed to explain the formation of 
disks and equatorial density enhancements around rapidly rotating stars. The idea of Wind 
Compressed Disks in nonmagnetic stars was developed by Bjorkman & Cassinelli (1993) and 
the concept of Wind Compressed Zones (WCZ) was presented by Ignace et al (1996). In 
the WCD model, trajectories of supersonic flow material from a rotating star travel to the 
equatorial plane. After the matter enters the disk region to form a WCD it continues to flow 
outwards. However, the angular speed of this material is far below the Keplerian orbital 
speed. An important similarity between the WCD model and the magnetically channeled 
model is that the fiow passes through a compression shock which leads to a significant 
enhancement of the density in the equatorial region. The WCZ model is a milder version 
of the WCD model and does not involve a magnetic field. Because it applies to stars that 
are rotating too slowly to form WCD structures, the WCZ model does not have a shock- 
compressed disk region and the enhancement of density is less than in the WCD model. The 
WCZ idea has been modified by Ignace et al (1998) to include the effects of a weak magnetic 
held. It is called the WCFields model. In this model, the magnetic flelds are unimportant 
in accelerating the flow but an enhancement of the fleld near the equatorial plane occurs 
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through compression of field lines. The WCFields model may be relevant in discussing the 
fiow in the equatorial region beyond and above the disk in magnetic rotator models. 

Babel & Montmerle (1997a) developed a magnetically confined wind-shock (MCWS) 
model to explain properties of Ap and Bp stars. In this model, a strong dipole magnetic field 
channels the fiow of wind from an early type star to form a disk. Babel & Montmerle (1997b) 

and Donati ct al (2001) have applied the MCWS model to hot stars in which observational 
data have been interpreted to indicate the presence of dipole magnetic fields with strengths 
in the range from 300 to 400 G and with significant obliquity. Unfortunately, these papers 
neglected the contribution of rapid rotation to the dynamics of the wind. The hot stars of 
interest to us, especially Be-type, are known to be extremely fast rotators and the effect of 
rotation on the dynamics of the wind will be more significant than that of the magnetic field. 
The MTD model addressed the question of disks in fast rotators with strong magnetic fields. 

Several phases in the development of disks around Be stars have been revealed by long 
term observational studies of these objects. Disks can form and disappear on time scales of 
about 30 years. When the disks of Be stars achieve a certain mass, as estimated by the equiv- 
alent width of the emission line, the stars tend to develop the observational phenomenon 
of Vj R variations. During the V/ R variation phase the star displays a modification in which 
the V component periodically changes so as to appear stronger and weaker relative to the R 
component. The V/R phenomenon has a typical time scale of about 7 years, which has been 
interpreted by Okazaki (1997) and others as the orbital period of a one-armed spiral pattern 
and is typically much longer than the orbital time of matter in the disk region around the 
star. 

In Section 2, we discuss the spatial features of flow regions or wind zones in our models. 
In Sections 3 and 4 we set up systems of equations that can be used to study properties 
of flows and winds in rotating stellar models with magnetic flelds. the models have axial 
symmetry and are in a steady state. In a general study of magnetic rotator winds it would 
be necessary to consider details of the magnetic field over the entire stellar surface. If the 
objective is restricted to studying the formation of disks, we need to consider magnetic fields 
that are sufficiently strong only over a smaller locahzed sector of the stellar surface through 
which fiux tubes emerge from the interior. 

In Section 5 we set up jump conditions that must be satisfied across the shock boundary 
between a fiow region and a disk region. We establish an important theorem in Section 6 
connecting the angular velocity of a disk region in which there is no meridional motion with 

the angular velocity at the photosphere. In Section 7 we consider conditions that must be 
satisfied during the fill-up stage in a disk region with no meridional motion and we derive a 
necessary condition for the formation of a Keplerian disk. In Section 8 we set up equations 
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that govern the structure of a disk region in a steady state. In Section 9 we discuss the 
rotation of disk regions during the fill-up stage and, also, the extent of Keplerian regions 
after a quasi-steady state has been reached. 

We consider two processes in the formation of Keplerian disk regions. These are a fill-up 
process and a process of adjustment of disk material with super-Keplerian rotation speeds 
into Keplerian orbits. We use the term quasi-steady Keplerian disk to describe a disk region 
that is formed after the processes of fill-up and adjustment into Keplerian orbits. We study 
models in which it is reasonable to consider the processes separately and we discuss this in 
Sections 9, 13 and 14. Several authors (Lee et al 1991; Porter 1999; Okazaki 2001) have 
suggested that viscous stresses play an important role in the structure and properties of 
disks. We assume that viscous forces are not important during the fill-up stage. 

In Section 10 we discuss some constraints on the magnetic field and in Section 11 we 
develop a prescription for a stream function and the corresponding dipole-type magnetic 
field configurations that enable us to study disk formation. We discuss further details of our 
models and develop a method of solution in Section 12. We present our results in Section 
13. Conclusions and possible interpretation of observations are discussed in Section 14. 

2. Flow Regions and Wind Zones 

We consider the flow of matter away from the stellar surface caused by traditional 
wind driving processes. While most of the flow from higher latitudes will travel to large 
distances away from the star in the form of a wind, material from lower latitudes can fiow 
into equatorial disk regions or even fiow back onto the star. In order to simplify terminology, 
we to refer to any fiow region surrounding the star other than a disk region as a wind zone. 
This includes a region in which matter fiows to a disk region or back towards the stellar 
surface. Figure 1 shows a schematic drawing of a meridional cross-section of the different 
wind zones and, also, the disk region during the fill-up stage. In each hemisphere of our 
models, we consider three wind zones. W^isk denotes the wind zone in which the matter 
fiows into the disk region. Wouter is the wind zone in which material actually leaves the 
stellar environment in the form of a wind. A third zone that may be present is Winner where 
matter fiows towards the equatorial plane and then back towards the stellar surface. When 
disks are present, we consider models with a single contiguous disk region D. 

We denote the photosphere by E^. Let Eq, which we call the initial surface, be a surface 
close to the the photosphere on which initial conditions for the fiow in the wind zones 
are specified. We can take Eq to be the sonic surface or some other appropriate surface. 
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Quantities associated with the photosphere and the initial surface are denoted by subscripts 
-k and 0, respectively. We identify different sectors of the photosphere and the initial surface 
that are linked to the different wind zones. At higher latitudes, we have the the sector 
Eo.ojtter of the initial surface from which material flows through the wind zone Wouter along 
open streamlines to a large distance from the star. The meridional speed of the flow in this 
zone reaches the meridional Alfven velocity at some distance from the star and exceeds this 
value at points beyond the Alfven distance so that the magnetic field hues are drawn out 
by the flow. The sector of the photosphere that is associated with So,outer is denoted by 
'^■k,outer- The flow from the sector Eo,disfc of the initial surface travels through the wind zone 
Waisk and passes through a shock surface T, shock to supply material to the disk region D. Let 
^*,disk be the sector of the photosphere that is linked to the sector So.disfc by the magnetic 
field. If a third wind zone Winner exists such that initially the fiow is towards the equatorial 
plane and then back towards the surface, wc take So^mner and S^^j„„er to be the sectors of 
the initial surface and photosphere that arc associated with Winner- If the sector S^^j„„er is 
present, it would be the closest to the equator. This zone does not play any significant role 
in our models. The meridional speed of the flow at all points in the regions Winner and Wdisk 
is less than the meridional Alfven speed so that the magnetic held is able to influence the 
flow in theses zones. Although we separate the flow into different zones, the basic equations 
that we set up in Sections 3 and 4 of this paper are valid for all flow regions or wind zones. 

An important feature in our models is the meridional streamline between the wind 
zones Wouter and Wdisk- We use the term limiting streamline to refer to this streamline, 
which separates the flow entering the disk region from the wind flow to large distances from 
the star. During the fill-up stage, the shock boundary, Ughock, between the zone Wdisk and 
the disk region D does not extend outwards beyond the limiting streamline. After the fill- 
up stage, the material in the disk region with supcr-Keplerian rotation speed diffuses out 
until a quasi-steady region is formed. The outer boundary of the quasi-steady disk region 
may extend beyond the limiting streamline and, if this happens, we assume that the wind 
material in the region Wouter ffows outside the extended disk region. 

We assume that, in addition to axial symmetry, streamlines of the flow and magnetic 
field lines have equatorial symmetry. Therefore, we need to set up our equations only in 
the Northern Hemisphere, which is determined by the direction of the angular momentum 
vector. Equatorial symmetry together with appropriate changes in sign, where necessary, 
will yield the values of quantities and ffow patterns in the Southern Hemisphere. In models 
with surface magnetic fields that have predominantly dipole-type structure we may assume, 
without loss of generality, that the field and the fiow in the Northern Hemisphere are in the 
same direction and in opposite directions in the Southern Hemisphere. 
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3. Equations for Flows and Winds 

In this section, we set up equations that are vahd for the material in all flow regions, or 
wind zones. We assume that the flow is steady, which is reasonable when all other phenomena 
of interest in the wind zones have time scales that are longer than the characteristic time 
of the flow. We consider a rotating hot star with a magnetic field. The star has mass M, 
luminosity L, equatorial radius it! and effective temperature Tg//. For simplicity, we assume 
that the star is approximately spherical and that the system is axially and equatorially 
symmetric. We use a stationary, inertial frame of reference with spherical polar coordinates 
(r, 6*, 4>). Let e^, eg, and be unit vectors in the r-, 9- and (p- directions. Subscripts m and 
denote meridional and azimuthal components of vectors. The velocity and magnetic field 
can be written in the form v = + v^e^ and B = + B^e^, respectively. Let Vm and Bm 
denote the magnitudes of Vm and Bm, respectively. The quantities Vm and B^ are always 
positive because they are the magnitudes of vectors. However, and B^ can be positive or 
negative because they are the 0-components of the vectors. 



3.1. Kinematic Conditions for Steady Flow 

The density p and the velocity v of the fiow material satisfy the continuity condition 

div(pv) = 0. (1) 

The divergence-free condition for the magnetic field is 

divB = (2) 

and the steady state MHD Induction Equation for material with extremely high electrical 
conductivity is 

curl(vxB) = 0. (3) 

Chandrasekhar (1956) and Mestel (1968) have obtained useful integrals of these differential 
equations. In an axially symmetric system, the azimuthal component of Equation (3) requires 
that 

v„ = , (4) 

where x = (vm-Bj„)/i?^ is a scalar function of position. Equation (4) shows that the 
meridional fleld lines coincide with the meridional streamlines. Also, Equations (3) and (4) 
give us 

= xB^+oursinO , (5) 
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where a; is a constant along a meridional streamline line. We interpret oj to be the constant 
angular velocity of the magnetic field line passing through the point Pq at the initial sur- 
face. Thus, if we picture the magnetic field line as being anchored at the point of the 
photosphere, then uj is the angular velocity of the star at P*. Equations (4) and (5) can be 
combined to give 

V = xB + ujr sin Q . (6) 

Equations (1), (2) and (4) require that the quantity defined by ,^ = must be constant 
along a meridional streamline in a hemisphere. This constant for a meridional streamline 
can be written 

i = PX = P ^5 ■ (7) 

Note that x > when and B„j arc in the same direction and x < when they are in 
opposite directions. In Section 4.2 we discuss properties of and P^ and show that the 
quantity XB4, is negative in both hemispheres. Therefore, Equation (5) implies that, relative 
to a stationary frame of reference, the rotational velocity of a fiow particle moving along 
a field fine is less than the co-rotation velocity a;r sin 6. If we consider points on a streamline 
through the point Po(ro, ^0, 0) on the initial surface. Equation (7) gives 

"r~ ~ ~R ' ^ 

where subscript refers to values at Pq. 

Equation (4) has important consequences for modeling stellar winds and disks. First, 
consider a star in which the dominant component of the surface magnetic field is an odd 
spherical harmonic, such as a dipole-type field. In both hemispheres, material fiows outwards 
from the stellar surface while the the magnetic field is directed outwards in one hemisphere 
and inwards in the other. Thus, x must have opposite signs in opposite hemispheres. This 
raises the question of continuity of field lines across the equatorial plane. We resolve this by 
noting that there are two possible scenarios for the flow of material close to the equatorial 
plane. In one of these, opposing flows from the different hemispheres will cross a shock 
surface and form a disk in which the radial velocity is zero or negligible. Here, shock 
conditions will account for the discontinuity in x- In this paper, we consider models where 
this condition prevails. In the other scenario with a dipole-type field at the surface, the flow 
near the equatorial plane is deflected parallel to the plane and wc do not have to consider 
the continuity of x across the plane. If the flow crosses a shock surface before moving 
parallel to the equatorial plane we expect a WCD region to be formed. In both scenarios, 
the flow near the equatorial plane in the region Wouter, outside the limiting streamline, will 
be almost parallel to the plane so that corresponding streamlines in opposite hemispheres 
are not connected. Part of this flow may form a WCZ. In the case of a star with a magnetic 
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field in which the dominant component at the initial surface is an even spherical harmonic, 
such as due to a quadrupole field, the field lines near the equatorial plane in the different 
hemispheres will be parallel to the plane and the question of the continuity of x does not 
arise. 



3.2. Dynamical Equations 

The general form of the steady state equation of motion for the flow of matter is, in the 
usual notation, 

/ ^ (l-r)GM curlBxB ^ 

p V • V V = - Vp - p ^ { e, + + Fune , 9 

where F = agL/ (iTiGMc) is the ratio of the stellar luminosity to the Eddington luminosity 
and (7e is the electron scattering opacity. The term involving F gives the force due to 
continuum radiation written as a fraction of local gravity. The force due to line radiation 
is Fiine — pgi^r: whcrc ql is givcu by the usual CAK expression (Castor, Abbott & Klein 
1975). We assume that the pressure is given by the isothermal equation of state 

P = Pcl, (10) 

which is a reasonable approximation for the wind zones and the disk. Here, is the isother- 
mal speed of sound. If we apply the continuity condition (1) to the left hand side of the 
^-component of Equation (9), and then multiply both sides by r sin^, we obtain 

A- f ■ n \ -a ^ curlBxB ^^ 

div^pvrsmfc^v,^ j = rsmt/ ^ ^ j . [11) 

The right hand side of this equation represents the torque exerted by the magnetic force on 
the flow material. This equation gives the integral 



rs\n9v^—y-^-^jrsui9B^ — C^ (12) 

where £ is a constant along a meridional streamline. We interpret C to be the total angular 
momentum per unit mass at points along the streamline. 

In general, the line radiation force, F/j^g, given by the CAK expression is not conser- 
vative. In models where this force is important, in order to obtain an energy equation we 
would have to flnd a suitably deflned potential function whose gradient will be approximately 
equal to Fune- Since our primary interest is in understanding the effects of magnetic flelds. 
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we consider models in which the hne radiation force is neghgible and we set Fune — 0. This 
is a vahd approximation in the case of late Be stars (Cassinelli et al 2002). We can write 
Equation (9) in the form 



12 2, (i-ncM 

-^^ + c^lnp-^ ^ 



, curlBxB 

pvx(Vxv)H . (13) 

47r 



We take the scalar product of this equation with v and use Equations (6) and (11) to obtain 
the BernouUi Equation 

^'2, (i-r)GM 

— + cjlnp-- ursme v^^ £ , (14) 

where £^ is a constant along a meridional streamline and has dimension of energy per unit 
mass. The term a;rsin^ expresses the work done by the magnetic torque on the flow 
material. If we move this term to the right hand side of the equation and compare the new 
equation with the corresponding equation for a nonmagnetic system, we see that the energy 
of the flow material is increased by the work done by the magnetic field. Only the azimuthal 
component of the magnetic force does work on the material and supplies energy to it. This 
is because the flow velocity, given by Equation (6), is the sum of a vector parallel to B and 
another in the azimuthal direction while the magnetic force is perpendicular to B. 

Our equations give us the set {^, u, C, S} of streamline constants for the flow of material 
in the wind zones. Here, ^ is the ratio of the meridional mass flux density to the meridional 
magnetic flux density and u is the angular velocity of a flux tube associated with a streamline. 
Also, jC is the total angular momentum per unit mass and £ represents an energy per unit 
mass and is not the total energy per unit mass of the flow material. 



4. Alfven Velocities and Azimuthal Components 

4.1. Alfven Velocities 

Since the magnetic field plays an important role in determining the nature of the fiow, 
it is useful to compare the fiow velocity with the Alfven velocity of the magnetic field, which 
is given in vector form by B/y^Trp. We introduce the quantities 

Am = —?= — and A<ty = —7= — • (15) 

Note that Am is always positive but Aife can be positive or negative depending on the sign 
of Bff,. Also, Am and \A((,\ are, respectively, the reciprocals of the meridional and azimuthal 
Alfven Mach numbers. 
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Streamlines flowing through the wind zone Wouter are associated with the material flow- 
ing away to large distances and are characterized by the presence of an Alfven point, Paif, at 
which Vm equals the meridional component of the Alfven velocity. Let subscript alf denote 
quantities at Paif- Then, we have 

4:Trpalf 

which satisfies the condition that Am,aif = 1- Equations (8) and (16) combine to yield 

Paif = 52 = ^2—' y^') 

-°m,0 "^m,0 

which is the density at the Alfven point on an open streamline and is determined by the 
conditions at the initial surface. 

In the wind zones Wdisk and Winners the magnetic field is strong enough such the con- 
dition Am > 1 is satisfied at all points in these zones. That is, the meridional component 
of the flow velocity is sub-Alfvenic in these zones. On the other hand, for a streamhne in 
the wind zone Wouter that emerges from the point Pq on So,outer) we have Am > 1 at all 
points between Pq and Paif, and Am < 1 for all points beyond Paif - The limiting streamline 
that separates the wind zones Wouter and Wdisk is the first meridional streamline above the 
equator that contains an Alfven point. Assuming that the system is equatorially symmetric, 
the Alfven point on the limiting streamline will be on the equatorial plane. 

For streamlines in all wind zones, we define two quantities Pc and which are constants 
for a given streamline but can be different for different streamlines. The first is given by the 
equation 

^c=72- = ^- (18) 

This may be interpreted as the characteristic density associated with a streamhne. The 
second quantity is defined by 

ujf^C. (19) 

We may interpret £ to be a characteristic distance from the rotation axis associated with a 
given streamline. On an open streamline in the region W outer i we have Pc — Paif and I is 
the distance of Paif from the rotation axis. For streamlines in the region Wdisk or Winner i 
we consider p^ and i to be streamline parameters that enable us to conveniently express 
properties of the flow. The value of £ for any streamline in the Northern Hemisphere is given 
by Equation (Al) in Appendix A. 
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4.2. Azimuthal Components 



If uj is the angular velocity at the point of the photosphere, we have f = u R sin6^. 
If 6^ = 6o, so that the points and Pq arc at the same latitude, Equation (5) gives the 
relation between the rotational velocities at the photosphere and the initial surface as 



1 - 



Equations (5), (7), (12) and (19) give us 

■l-(Pc/p){^V(^'sin2^)} 



(20) 



and 



— Lor sin 9 



— 47rpx cur sin 9 



1 - (Pc/P) 

l-{ey{r^sm^9)} 



(21) 



(22) 



1 - (Pc/p) 

Because p > Pc when r sin ^ < £ and vice versa, Equation (22) requires that and x must 
have opposite signs when a; > 0. That is, is always negative. In our models, x is 
positive in the Northern Hemisphere. Hence, B^ and Ac/, are negative in that hemisphere. 
Expressions for azimuthal components and Alfven velocities that are useful in applications 
to models are given in Appendix A. 

The equations that we have set up are valid for any appropriate rotation law at the 
stellar surface. We introduce the rotation rate parameter a^, — ai,{9^) for the star at the 
point at the stellar photosphere by taking uj — cucrit where 



^crit — 



GM{1-T) 
R^ 



1/2 



Then, the rotational speed at the point of the photosphere is given by 

v<p,* = u}critCt*Rsin9i, . 
Prom Equations (20) and (24) we obtain 



V(l),0 — i^crit 



ro sin 9q 



(23) 



(24) 



(25) 



for the rotational velocity of the flow material at the point Pq of the initial surface. In models 
where the angular velocity is constant over the sector T,-k^disk of the photosphere, we write 
= a, where a is a constant. 
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5. Jump Conditions across the Shock Boundary 

The material from the sector So,disfc of the initial surface passes through the shock 
surface T^shock before entering the disk region D. We consider jump conditions across T^shock 
during the fill-up stage of the disk. We assume that the shock is isothermal, which is a 
reasonable approximation for the types of models that we consider (CassineUi et al 2002). 
In models where the dominant component of the stellar magnetic field is an odd spherical 
harmonic, such as a dipole-type field, we assume that there is a region of the disk that has 
no meridional motion during the fill-up stage. We refer to this pure rotational region of 
the disk as a pre-Keplerian region and we expect that this region or part of it evolves into 
a quasi-steady Keplerian disk region after the fill-up stage. On the other hand, in models 
where the dominant component of the stellar field is an even spherical harmonic, we expect 
that the material in the disk region will have an outward fiow along field lines that are 
parallel to the equatorial plane. For convenience, we refer to such a region as a WCD region. 
Another possible situation is one in which matter in the inner part of the disk region is 
pre-Keplerian and the outer part is of WCD-type. Although our main interest in this paper 
is in pre-Keplerian regions, since one of our objectives is to set up equations with general 
validity, we include conditions for Pre-Keplerian and WCD regions. 

In each hemisphere, T^shock extends from the inner boundary of the wind zone Wdisk to 
its outer boundary. In the northern hemisphere of the meridional plane (f) — 0, the inner 
boundary of Wdisk coincides with the streamline from the point Po,mn on the initial surface 
to the point Qinn on T,shock- Let Xj„„ be the point on the equatorial plane that is associated 
with Qinn- The outer boundary of W^isk is the limiting streamline from Po,/im on the initial 
surface to the point Qum on T^shock and is associated with the point Xu^ on the equatorial 
plane. Let Po{ro, 9o,0) be a general point on the sector So,disfe so that 6'o,/im < 9o < Oo^inn- 
The meridional streamline from Pq travels to the point Q on the shock surface, which is 
associated with the point X(rx,7r/2,0) on the equatorial plane. In models where the disk 
region is a combination of a pre-Keplerian region and a WCD region, the inner part of D 
from Xinn to the point X^jcd is pre-Keplerian and the outer part from Xy^cd to is a 
WCD region. The point Xwcd{rx,wcd, 7r/2, 0) that separates the pre-Keplerian and WCD 
parts of D is associated with the points Po,wcd{ro, 6*0 .u,cd, 0) on T^o^isk and Q.ujcd on T.,hock- 
The equations that we set up in Section 3 for the wind zones are also valid in a WCD region, 
where there is a steady radial fiow of material superimposed on the rotational motion. 

Let n be the unit normal vector to T^shock from the wind zone into the disk. We define 
the unit tangent vector t in the meridional plane as t = nXe<^. We use subscript n to denote 
components of vectors in the direction of n and the notation A[...] to represent the change 
in a quantity across the shock boundary We use a hat to denote quantities in the post-shock 



disk region and the subscript fill to indicate the values of quantities at the end of the fill-up 
stage. 

In a region of the disk with no meridional motion, the mass increases steadily. In such 
a region, the density and the height will change with time until a steady state is reached. 
The magnetic field is assumed to be steady during the fill-up stage of the disk. Let h be the 
height of the shock surface at a radial distance r along the equatorial plane. Let 5 be the 
angle between the tangent line to the meridional cross-section of T^shock and the equatorial 
plane. In general, 6 changes with r and with time, and we assume that it is a small positive 
angle. Then, h increases with r so that dh/dr = tan 5 > 0. Although h and d change with 
time, they do so on a time scale of the order of the fill- up time tfm of the disk, which is 
much longer than the fiow time tfiow = R/vm- 

In both the pre-Keplerian and WCD regions of the disk, conservation of magnetic fiux 
requires that the normal component of the magnetic field be continuous. Thus, 

A[S„]=0. (26) 

Also, the tangential component of the electric field should be continuous across Y!i shock- The 
electric field can be written as E = — vxB/c. Therefore, the jump condition requires that 
the tangential component of vxB should be continuous across the shock surface. The 
tangential component of vxB has contributions from its ^-component and its t-component. 
The 0-component of vxB is equal to VmXBm, which is equal to in a flow region and 
a WCD region because of Equation (4), and in a pre-Keplerian region because = 0. 
Therefore, for both pre-Keplerian and WCD regions the jump condition reads 

A[(vxB) -t] = 0. (27) 



The boundary conditions for mass flux are different for pre-Keplerian and WCD regions. 
We consider them separately. In a pre-Keplerian region, there is no radial motion and 
material flowing into such a region will steadily increase the mass in that region. This gives 

P^;„sec5=^. (28) 

Because 5 is quite small in thin disks and we can take sec 5 ^ 1. Also, as the disk region 
continues to flll up, pq becomes much larger than pq and the time taken by h to change 
signiflcantly becomes quite long compared with the time scale t^iow for the flow of material 
towards the disk. 

The continuity of mass across the shock surface into a WCD region gives 



(29) 
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Then, for points X between X^^cd and Xum such that ^o,Hm < < ^o,wcd, since any inflow 
from the pre-Keplerian part of the disk into the WCD part is neghgible over the time scales 
of interest to us, mass conservation for flow into a WCD region requires that (Bjorkman & 
CassineUi 1993) 

\,wcd 

2nrQ sin 9 po 'Vrfi d9 — 2nh rx pVr , (30) 



where Vr is the speed of the WCD flow parallel to the equatorial plane averaged over the 
height of the disk. In models where po and are constant, this condition simplifles to 

Po Vr,o (cos 00 - cos 9o,y,cd) ^hrx pvv (31) 

Equations (30) and (31) give the rate of mass flow in the Northern Hemisphere only. The 
total rate of mass flow is twice the value of the integral. 

If we substitute curlBxB = (B ■ V)B - V(i?^)/2 in Equation (9), the continuity 
condition for momentum flux across T^shock into a pre-Keplerian or WCD region yields 



S„B 



n V + \ p+ — n - 



0, (32) 



Stt / 47r 

where we have assumed that the line radiation force is negligible. 

Next, we consider the conditions that must be satisfled by the streamline constants i^, a;, 
C and £ for flow from the wind zone into a WCD region across the surface shock- Equations 
(26), (29) and (7) imply that 

A[e] = 0. (33) 
Also, Equations (6), (27) and (26) together with the fact that e<^Xt = n require that 

A[u;] = 0. (34) 

Equation (12) and the azimuthal component of (32) ensure that 

A[£] = 0. (35) 



In general, the streamline constant £ in the Bernoulli Equation is not continuous across 
the shock boundary of the disk region. We do not consider jump conditions for streamline 

constants across the shock boundary of a pre-Keplerian region because there is no meridional 
motion there. However, since oo is associated with the rotation of a magnetic field line, in 
Section 6 we discuss the continuity of u) across the shock boundary of a pre-Keplerian region. 



- 15 - 



6. A Theorem on Disk Rotation 

Here, we consider the jump condition given by Equation (27) at a point Q{rQ,9Q,0) 
on the shock boundary of a disk region with no meridional motion. In the pre-shock the 
flow region, Equation (6) gives (vxB) ■ t = —urq sin 9q Bn^g ■ In the disk region we have 
V = v^e,^ and this gives (vxB)-t = —'V(f,,Q Bn,Q. Also, Equation (26) states that Bn,Q — Bn,Q. 
Therefore, Equation (27) gives 

Vcj>,Q ^ujrq sin 9q . (36) 

This equation shows that the point in the disk region that is adjacent to Q rotates with 
angular velocity cu. In Section 3 we noted that the angular velocity cu associated with a field 
hne is the same as that of the point P^, of the star where the field hne anchored. Thus, we 
have an extremely important result for the rotation in disk regions where the meridional 
speed is negligible. We state this in the form of a theorem as follows: 

Theorem 1 Let Q be any point on the shock boundary between a magnetic rotator flow 
region and a disk region of an axially symmetric star that is in a steady state. Suppose 
that there is no meridional motion in the disk region and that the normal component of the 
magnetic field in the disk region at Q is not zero. Let Q be the point of the disk region 
adjacent to Q on the field line through Q. Then, the angular velocity of the disk region at Q 
is equal to the angular velocity, cu, at the point of the star at which the magnetic field line 
through Q is anchored. 

Furthermore, if the system is equatorially symmetric and the magnetic field line through Q is 
continuous across the disk, then the angular velocity at all points of the disk region along the 
field line will be the same as that at for time periods over which the disk region remains 
steady. 

The first part of the theorem follows from Equation (36) and requires that the meridional 
speed in the disk region be zero. To establish the the second part we use the equation for 
the pre-Keplerian region that corresponds to Equation (5). Since = 0, we have x = 
and we obtain i)^ — cjrsind, where a) is constant along a meridional field line. Comparing 
this equation for t)^ with Equation (36) we see that u = cu for the field line through Q. Thus, 
all points of the pre-Keplerian region along a field line will have the same angular velocity 
as the point of the star where the field line is anchored over time scales where the region is 
steady. During the fill- up stage, the fiow time, tfiow, is much shorter than the the fill up time, 
tfiii, so that over time scales of the order of tfiow the second part will be true. A situation 
in which the second part of the theorem will be true over longer periods of time that are 
comparable to tfm is when the region 'Ei.^disk is in uniform rotation. Then the corresponding 
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pre-Keplerian region of the disk will have the same constant angular velocity as that sector 
of the star. The models that we consider in our applications satisfy this condition. 

We emphasize that our theorem is quite different from Ferraro's Law of Isorotation 
although part of it may appear to be similar. We require the presence of a shock surface and 
in the region between the stellar surface and the shock surface there is meridional motion as 
well as differential differential rotation of the outflow material. In models satisfying Ferraro's 
Law, we can easily show that there should be no meridional motion outside the stellar surface. 
From our discussion in Section 3, wc know that x must be constant along field lines and 
have opposite signs in opposite hemispheres. If there is no shock surface then the continuity 
of X across the equatorial plane requires that x = along every field fine. Therefore, 
must be zero along every field line and there can be no meridional motion anywhere in an 
atmosphere where Ferraro's law is satisfied. 

7. Conditions for the Formation of Pre-Keplerian Disk Regions 

In applications to models, wc consider conical shock surfaces given by = (7r/2) — 5 in 
the Northern Hemisphere, where 6 is small and has the same value for all r. It may increase 
steadily with time during the fill-up stage. In such a model, the cross-section of T^shock is a 
straight line given hy h — r tan 5 and we have n = e5t and t = e^ . The jump conditions that 
we set up in Section 5 give us the following equations at the point Q on the shock boundary 
of a pre-Keplerian region of the disk, where Vr — Q and vq — 0: 




(37) 





(38) 









1 



(40) 



(39) 



Equation (40) can be combined with Equation (12) to obtain 




(41) 



At the end of the fill-up stage, if pfm is the disk density averaged over the height of the 
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disk at any point, we can use the jump condition (32) for momentum flux to write 

•c2 + y2\ (B'-B 



Pfill,Q 



(42) 
Q 



Equation (39) shows that B^q is shghtly less than B^q and Equation (40) shows that q is 
slightly larger than B'^ q. In stellar models with dipole-type flelds, we have B^^q ~ near the 
equatorial plane and because of Equation (38) we can write B^ — B^ — B'^ q — B'^^q, which 
will be negative. Hence, in such models, the density of a Keplerian disk will be somewhat 
less than that in a nonmagnetic model. Also, if the magnetic fleld present is too strong, 
Keplerian disk formation may not be possible. A necessary condition for disk formation can 
be obtained by requiring that pfm > in Equation (42). Combining this with Equation 
(40), using the relation i?| = Anpv^A^ and noting that Vq ^ Cs near the shock surface, we 
obtain a condition for Keplerian disk formation as 



V0,Q ^ ( 1 - 2Am,QA^^Q 



1/2 



(43) 



That is, for a given rotational velocity of the pre-shock flow at a point of the shock surface, the 
normal component of the flow velocity at that point must be larger than a value determined 
by the magnetic fleld. Also, if and v^^q are directly related to the wind terminal velocity, 
Vooi and the angular velocity, a;, of the star , respectively, we expect that for a given value 
of Voo there will be an upper bound on uj beyond which a Keplerian disk will not form. 

We deflne the flU-up time at radial distance r for a pre-Keplerian disk region, to be 

hfiii pfiii , . 

Our expression for tfm is different from the one used by Cassinelli et al (2002). Since 5 is 
small during the flU-up stage for thin disks, we take cos 5 1. Then, Equations (28) and 
(44) give us tfm — hfm pfm/ pvn- For our conical shock surface with hfm — r tan Sfm, we 
have 

fpfm\ f r tan 6 fm \ . . 

where the value of pfm is given by Equation (42). After the fill-up stage, the density of a 
Keplerian disk region will be sufficiently large and we assume that the disk temperature will 
be determined essentially by stellar radiative heating. 
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8. Equations for Disk Regions 



In this section we consider equations for a disk region D that is in a quasi-steady state. 
We assume that f e = in D and the equations for the disk region are written for values of 
quantities averaged over the height of the disk. The r- and 6- components, respectively, of 
the momentum equation for the disk are 



and 



Here, 



dVr 

dr 



Vr d 

7 Yr ^^"^^ 
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p dr 
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r sin 9 dO 



— (r smOB^) 



(46) 

(47) 

(48) 
(49) 

(50) 



The last term in Equation (47) represents the effect of viscosity (e.g., Okazaki 2001) and v 
denotes the viscosity coefficient, which has a value of order 0.1. The forces per unit mass 
due to the magnetic field in the radial and azimuthal directions are given by !Fr and !F^, 
respectively. The quantity fr is the resultant of the three main forces (per unit mass) that 
control the radial momentum of the disk. These equations are supplemented by a mass 
conservation equation that has the form AnrhpVr = Mdisk, which represents the amount of 
mass flowing across a cross-section of the disk at radius r. In WCD regions, this equation 
can be replaced by Equation (30). In a strict Keplerian disk Vr and Mdisk are zero. 

In a Keplerian region of D we assume that Vr and Vp are negligible. In stars where the 
force due to line radiation is not significant we take Fune = 0. Then, Equation (46) requires 
that = and Equation (48) gives us an expression for the Keplerian rotation speed, ■0</,,k, 
in the form 



(51) 



where ( — r'^Tr/ {GM[1 — F)} is the ratio of the radial component of the magnetic force to 
effective gravity. Note that C, becomes negative when Tj. is directed inwards and C ~ when 
the magnetic force in the disk is small compared to effective gravity. Also, in a Keplerian 
disk region Equation (47) requires that Tcf, — {d (r'^uhpcg) /dr}/{hpr^) — 0. The results 
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presented for our models in Section 13 show that the magnetic force and the viscous force 
become neghgible at the end of the fill-up stage. Thus, both terms on the right hand 
side of Equation (47) become negligible. Although we do not study the evolution of the 
disk region after a quasi-steady Keplerian region is established, if a viscous force becomes 
significant, then in order to maintain Keplerian rotation, either the magnetic torque on the 
disk material should offset any decrease in angular momentum due to viscosity or there 
should be a continuous supply of angular momentum to the disk material. In the case of a 
WCD region, the steady state equations are essentially the same as those that we derived in 
Sections 3 and 4 for the wind zones with the additional conditions that vg — and Bg — 0. 



9. Rotation and Extent of Disk Regions 

We introduce the non-dimensional variable x = rx/R to represent the radial distance 
of the general point X{rx: 7r/2, 0) on the equatorial plane. We consider a pre-Keplerian disk 
region with shock surface given by ^ = (7r/2) — 5 in the Northern Hemisphere. Suppose 
that Q{rQ, {n/2) — S, 0) is the point on the shock surface such that the meridional field line 
through Q passes through X. When the disk is thin we assume that Q is vertically above 
X. Then, we have rx = tq cos 6. This is a good approximation in the case of a thin disk 
with a dipolc-typc field because the field line through X is perpendicular to the equatorial 
plane at that point. Wc assume that all points in a narrow vertical column through X in 
the disk region have the same angular velocity. 



9.1. Rotation of a Pre-Keplerian Region During the Fill-up Stage 

Theorem 1, which we derive in Section 6, gives 

v^^ur^ LOcrit a* r (52) 

for the rotational velocity of disk at distance r from the center of the star during the fill-up 
stage. The angular velocity, u, is constant along a field line and can have different values 
along different field lines. Thus, along the equatorial plane, cj is a function of r. To obtain 
this function wc must know c<j as a function of colatitude 6 at the photosphere and, also, 
the stream function, t/'(r, 6'), for the field lines linking the points P^, on the sector T,^^disk to 
the points X in the disk region D. However, even when the stream function is not known, 
we are able to discuss the rotation of a pre-Keplerian disk region of a star where the sector 
'^*,disk has an approximately constant angular velocity cu. In this case, the corresponding 
pre-Keplerian disk region will have the same constant angular velocity. Before we consider 
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models with simple rotation laws at the photosphere, we consider the question of whether it 
possible for exact Keplerian rotation to be established in a pre-Keplerian region at the end 
of the fill-up stage. We look at two examples which show that this is very unlikely in real 
stars. 

In a general model, the condition for Keplerian rotation of the disk is obtained using 
Equations (51), (52) and (23). It reads 

x^al^l-C, (53) 

where C is function of r and a* is a function of 6'*. Suppose that we consider an idealized model 
in which the angular velocity distribution at points in the sector Tii^^disk of the photosphere 
results in a pre-Keplerian disk region that has exact Keplerian rotation. Although such a 
possibility may be mainly of theoretical interest, it is quite instructive. Also, suppose that 
the disk is at the end of the fill-up stage when C, has become negligible. In such a model. 
Equation (53) requires that = x~^^^ for exact Keplerian rotation at a point X of the disk. 
Suppose that all points of the disk from Xi to X2 are in exact Keplerian rotation and let P^ i 
and P*^2, respectively, be the points at the photosphere that are linked to Xi and X2 by the 
magnetic field. Let us select xi = 1.3 and X2 = 3.5. Then, the rotation rate at points in 
the sector E^^^isfe of the photosphere must decrease from 0.67 at P^^i and 0.15 at P^^2- We do 
not expect that the rotation rate at the photosphere in real stars will decrease so drastically 
for a relatively small increase in latitude. Thus, for realistic angular velocity distributions 
at the photosphere, we expect that if the rotational velocity of the pre-Keplerian disk has 
a Keplerian value at a point X^^p then the rotation rate will be super-Keplerian at points 
between Xkep and Xu^- 

Another hypothetical model with exact Keplerian rotation is one in which the magnetic 
field in the disk region is quite strong and has a configuration such that the magnetic force 
is able to offset the difference between the centrifugal force and gravity. However, for the 
types of models that are of interest to us, the magnetic field required will be far too strong. 
This is because the radial component of the magnetic force in the disk must satisfy the 
condition ( = 1 — x^ a"^ and, for example, if a star with a = 0.6 is to have a pre-Keplerian 
disk extending to a distance oi x — 3.5, the equations that are given in Section 8 together 
with the disk density at the end of the fill-up stage given by Equation (42) show that the 
value of C must be approximately equal to —14. Such a large value of C in the disk region 
will require an exceptionally strong magnetic field at the photosphere. Also, the topology of 
the field would have to be highly contrived to produce a magnetic force in the disk region 
that has the required magnitude and direction. 

The different cases that we have considered here suggest quite strongly that the most 
plausible scenario for the formation of a quasi-steady Keplerian disk region is one where a 
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pre-Keplerian disk region is formed with material having super-Keplerian rotation speeds 
and then this material diffuses into Keplerian orbits when the magnetic force becomes small. 



9.2. Radial Extent of Quasi-steady Keplerian Disks 

Consider a model in which the sector T^i^^disk the photosphere is in uniform rotation 
with angular velocity cu. Let the surface magnetic field be sufficiently strong for the formation 
of a pre-Keplerian region of the disk during the fill-up stage. If there is no meridional motion 
in the disk region, our theorem implies that the angular velocity of this prc-Kcplcrian disk 
region will be constant and equal to u. Suppose that the surface field is not too strong and 
that at the end of the fill-up stage the disk density pfm becomes sufficiently large such that 
the magnetic force in the disk becomes negligible compared with the centrifugal force. Let 
Xkep be the innermost point of the disk where the rotation speed is Keplerian. To simplify 
our discussion, we take X^^p to be at the inner boundary of the pre-Keplerian region of the 
disk. Let Xkep be the value of x at the point Xkep- 

If a is the rotation rate of the sector T,^^disk, Equations (24) and (53) give 

Xkep = a"^/^ . (54) 

The rotational speed at points in the pre-Keplerian region between X^ep and Xum will exceed 
the local Keplerian speed. Let X^g^ be the point at which the rotational speed of the disk 
reaches the escape velocity so that Xesc = 2^/^q;~^/^. In some models, Xgsc will be located 
between X^ep and X^^ and in other models the rotational speed may not reach the escape 
speed between Xkep and Xum. Let Xint denote the nearer of Xesc or Xum to the point Xkep- 
Since the magnetic torque on the disk material is negligible after the fill-up stage, we expect 
that the super-Keplerian material in the region from X^ep to X^nt will diffuse outwards into 
Keplerian orbits while conserving angular momentum if viscosity is negligible. The material 
in the region beyond X^nt will have sufficient kinetic energy to flow away from the star. 

When a quasi-steady Keplerian disk region is formed after diffusion of the material in 
the pre-Keplerian region, let the point Xend be at its outer boundary. Then assuming that 
the material at X^nd arrived from Xi^t through angular momentum conservation, we have 

^end ~ •^int ^ • ThuS, 



24/3^-2/3 iixu^>X 



^end — S 4 2 -f ^ V"^ / 

In order to form a disk region with Keplerian rotation we require that > x^ep- When 
Xiim is slightly larger than Xkep the Keplerian region will be in the form of a ring. For larger 
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values of xum the radial extent of the Keplerian region increases until xum — Xesc, which 
we refer to as our optimal model When xum increases beyond Xesc, there is no increase 
in size of the Keplerian region. The radial extent of the largest Keplerian disk region is 
given by < a; < 2^/^ a^^/"^. If viscosity plays a role during the process of diffusion of 

super-Keplerian material into Keplerian orbits after fill-up, the radial extent of the Keplerian 
region will be slightly different. The presence of a weak viscous force may be helpful in the 
process of readjustment of angular velocity. 

We should point out that there are other situations in which quasi-steady Keplerian 
disks with different radial extents can be formed. For example, if the magnetic fiux tube 
that assists in forming the disk region is localized over the sector T^i^^disk such that the rotation 
speed at the innermost point X^^n of the pre-Keplerian disk region is super-Keplerian, then 
the the inner radius of the quasi-steady disk that is formed will be larger than oT'^I^R. 



In this section we derive expressions for Am,o,iim ^-nd .4<^,o,Hm that are useful for finding 
solutions. In our models we require that Am,o > 1 at the points Pq on the initial surface. In 

the wind zone Wdisk, the streamline from Pq travels to the point Q on the shock surface and 
Am > 1 at all points between Pq and Q. In the wind zone W outer: the streamline from Pq 
travels into space and Am > 1 at all points between Pq and Paif-, which is the Alven point on 
the streamline. In Section 4.2, we found that Bff, <{) and A(f, < in the Northern Hemisphere 
in our models. Equations (A4) and (A5) show that A^, will be negative if Am > (i/{rsin6) 
because £^ > sin^ d for points between Pq and Q, or between Pq and Paif. Therefore, the 
field at the initial surface should satisfy the condition Am,o > £/{rosm9o), which requires 
that Am,0,lim > Viim, where 



The dynamical equations governing the flow will also impose constraints on the flclds at the 
initial surface but we do not consider them here. When the density distribution in the wind 
zone is known. Equation (18) gives 



10. Constraints on the Magnetic field 



10.1. Magnetic Field at the Initial Surface 



yiim 



Xq sin OQ^um ' 



(56) 




(57) 



This is not an explicit equation for Am,o,iim because the value of px,iim depends on the 
value of Am,o,iim- Equation (12) gives a steady state solution of the azimuthal equation of 
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motion, which leads to the expression for Acj) in Equation (A5). When an appropriate value 
of Am,o,iim has been determined, the value of A^^q^u^^ is given by 

Acj>,Q,lim = Am,0,lim i ^^"^ 2 ) • (^^) 

\ m,0,Um yiim J 

Equations (57) and (58) are constraints that must be satisfied in a magnetic rotator wind 
model in a steady state. 



10.2. Critical Values of the Surface Magnetic Field 

Our approach to Keplerian disk formation does not depend on which theory of merid- 
ional circulation speeds is valid and we consider circulation only as a possible cause of some 
of the time-dependent behavior of disks. The radial and azimuthal components of the surface 
magnetic fields in rapidly rotating hot stars are subject to certain bounds (Maheswaran & 
CassineUi 1988, 1992), which are 



and 

"47rG'M(l-r)(l-Q;2)p^ 



1/2 

(60) 



R 

Subscript circ denotes quantities relating to meridional circulation. The effect of rotationally 
driven circulation on meridional magnetic fields have been studied by Maheswaran (1969) 
and by Roberts & Wood (1985). The critical value given in Equation (59) is important for 
our models. If the radial component 5^,* of the surface magnetic field is less than -6^,*,^, 
then meridional circulation currents driven by rapid rotation can draw the meridional field 
lines beneath the stellar surface in a time of order tdrc — '^T^fdrc/vdrc- In this case, the 
meridional magnetic field that provides the torque on the fiow material will not persist over 
periods of time that are much longer than tdrc- However, if tfm is less than tdrc a disk can 
form before the field lines submerge. For circulation to be involved in the time variation of 
disks, the meridional fields at the surface should be strong enough to assist in the formation 
of disks but weaker than -Br.,,t,L- Also, a meridional component of the field should be able to 
periodically emerge from the stellar surface after being submerged by the circulation. The 
fields need not cover the entire photosphere but emerge at least through the sector S^^^isfc- 
The causes for the resurgence of meridional field lines could be the buoyancy of magnetic 
fiux tubes or other MHD phenomena that occur in the stellar interior. We know that the 
Sun displays such phenomena and it is reasonable to expect that other stars do the same. 
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The critical value given by Equation (60) is used for ensuring that B^^^ is less than B^^^^u 
so that the sum of the centrifugal force and magnetic force does not exceed gravity at the 
photosphere. However, the magnetic fields in our models easily satisfy this condition. 

Meridional circulation velocities in rotating stars have been discussed by several authors, 
e.g.. Sweet (1950), Baker & Kippenhahn (1959), Maheswaran (1968), Pavlov & Yakovlev 
(1978), Tassoul & Tassoul (1995). Unfortunately, no circulation models are available for 
rapidly rotating hot stars with winds. Sweet's first order expansion method gave relatively 
modest values for circulation speeds near the stellar surface. Maheswaran showed that 
the first order expansion method was inadequate when the angular velocity is taken to be 
constant because the term in 1/p does not appear in the first order term for the circulation 
speed. In the case of uniform rotation, the 1/ p term appears in the second order and higher 
order terms so that the circulation speeds are much larger than the values estimated by 
Sweet. For nonuniform angular velocity distributions, the 1/ p factor in the circulation speed 
appears in the first order term (Baker & Kippenhahn 1959; Pavlov & Yakovlev 1978). Tassoul 
& Tassoul (1995) find that when eddy viscosity due to turbulence is included, the 1/p factor 
disappears and circulation speeds are similar to those obtained by Sweet (1950). However, 
the Tassoul & Tassoul model is not appropriate for hot stars with winds because it imposes 
the boundary condition that v.n = 0, which is essentially the same as enforcing Vr^circ = 
at the stellar surface. The fact that they obtain values of the circulation speed that are 
smaller than the velocities for mass-loss rates near the surface implies that their boundary 
condition is not appropriate for hot stars with winds. Also, they have not treated the surface 
region in detail so that it is not clear whether the choking of circulation by eddy viscosity 
is a consequence of the stringent boundary condition which requires the fiow to become 
transverse at the surface. The boundary condition used by Maheswaran (1968) involves only 
the conservation of mass so that Vj.^circ does not have to be zero near the photosphere. In 
the absence of models with more appropriate boundary conditions, we use the circulation 
velocity that is consistent with the results of Maheswaran (1968) and Pavlov & Yakovlev 
(1978), which is 



Here, fl is the angular velocity at the equator and Ail is the change in angular velocity from 
equator to pole, p = 3M / (AttR^) is the mean density of the star. In the case of stars that 
are in almost uniform rotation we can write the circulation velocity in the outer part of the 
radiative diffusion zone in the form 




(61) 



V, 



arc 



3L(l-r)^a^ 
AnpdrcGMR' 



(62) 
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We should point out that even if the Tassoul & Tassoul approach to meridional circulation 

is appropriate, the theory that we propose for Keplerian disk formation is not affected. The 
only difference is that meridional components of surface fields with the strengths that we 
require will be able to survive for much longer periods of time because the circulation speeds 
are very small and the turnover time is long. In this case, rotationally driven circulation will 
not be a factor in the time variation of disks. 



10.3. Time Scale for Ohmic Diffusion in the Disk 

The azimuthal component 5^ of the magnetic field in the disk region will be in opposite 
directions in opposite hemispheres and we wish to determine whether there could be early 
decay of the field through Ohmic diffusion because the height of the disk is small. The 
diffusion time scale of the field is given by tdi/f = ^TrL^/r], where L is the length scale of 
the field and rj is the magnetic diffusivity. Using the expression given by Spitzer (1962) for 
diffusivity transverse to the field, we obtain rj — 1.42 x lO^^T""^/^ cm^ s~^ for disk regions. 
Then, Equation (45) for the disk fill-up time together with L = hfm, gives 

tfiii ^ 1.13 X 10^^ pfiu ^gg^ 
tdiff r pvn T3/2 tan 6 fm 

We find that the diffusion times for the fields in the disk regions of our models arc very much 
longer than the fill-up times so that decay due to Ohmic diffusion will not be a factor. 



11. Dipole-type Magnetic Fields in Wind Zones 

Since we wish to study the formation and properties of Keplerian disks, we consider 
models in which meridional fields have a dipole-type structure. In particular, we look for 
solutions with stream functions that have the form given by 

^, sin^ ^ sin^ ^0 / ,n 

= (64) 

where A is a constant. The field lines that coincide with streamlines given by this equation 
correspond to a dipole-type field with 

2C AC 
Br = „ , , cos 6 and Bq ~ sin 9 , (65) 

where C is a constant. When A = 1 this gives a regular dipole field. However, even if the 
magnetic field at the stellar surface is that of a regular dipole, the field lines in the wind 
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zones may be slightly modified by the flow. In realistic models it is likely that the values of 
A will be smaller than 1. When A approaches the transverse component Bg approaches 
and when A = the field becomes purely radial. We consider models in which < A < 1. 

In our applications, the disk has a shock boundary given by ^ = (7r/2) — 5 in the 
Northern Hemisphere. If the meridional streamline from the point Po{ro, do, 0) on the initial 
surface arrives at the point Q{rQ, (7r/2) — S, 0) on the shock boundary, we have 

r, = rj!^Y\ (66) 

As described in Section 9, we take X{rx, '7r/2, 0) to be the point on the equatorial plane that 
is associated with the point Q such that rx = tq cos 5. Equation (66) gives the relationship 
between 9q and x in the form 

sin2^o= C^)\os^^H. (67) 

The initial surface in our models is close enough to the stellar surface so that we can take 
ro ~ R. Also, the disk is sufficiently thin such that cos 5 ~ 1. Then, Equation (67) can be 
written as ^ 

«i^^o = P^- (68) 

We can construct models by prescribing the values of any two of the quantities xum-, 0o,iim. 
and A. Then, Equation (68) fixes the value of the third quantity. In models where Oo,Hm and 
xiim are known, we obtain 

^^_21n(singo,,J 

Equations (65) and (67) yield 

f A2 ^ 

Bm,Q = S^,0 I [4 (^A _ 1) + X^] / " ^^^^ 

We assume that over the sector "Eo^disk of the initial surface, the density, po, and meridional 
speed of the flow, Vm,o, are constant. Then, using B^fi — V^'^Po Vm,Q Am,o we obtain 

1/2 



^Hm [4 jx' - 1) + A^] 

[4 (xL - 1) + 



(71) 



For our dipole-type fields, we obtain yum — for the quantity defined in Equation 

(56). This expression can be substituted in Equation (58) to compute the value of A^fi^Hm 
in terms of xum and Am,o,iim- 
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12. Models and Applications 

If we know the density, velocity and magnetic field strength at the stellar surface, we 
can use our system of equations and boundary conditions to compute solutions for the wind 
zones and the disk region. Because the flow streamlines and the magnetic field lines coincide 
in a meridional plane, a convenient approach is to introduce a meridional stream function 
■0(r, 9) for the wind zone and obtain solutions for it in terms of known values at the initial 
surface. In this paper we do not attempt to completely solve our system of equations. For 
the present, we focus on the formation and properties of Keplerian disks. We derive two 
results that help us develop a method to pursue this goal. For given values of 6q and ^^,0; 
Equation (Al) gives di/dAm,o > 0, which means that when -Bm,o increases so does xnm- 
Also, for given values of 6^0 and Am,Q, we obtain di/dAcpfi < 0. Since A(j,,Q < 0, when |i^</,,o| 
increases xum also increases. 

We consider models in which ijj{r. 9) has the form given by Equation (64) and compute 
solutions for a large number of different cases. First, we select a value for the uniform 
rotation rate, a, of the sector T,^^disk- Then, we select a value for ^^o,«jm from a range of 
plausible values. This specifies the location of the point Po^um, which is the upper boundary 
of the sector So,djsA; on the initial surface. Next, we choose an appropriate value of xum, 
which fixes the magnitudes of Bm,o,Hm a-nd S<^,o,iim at Po,iim- In our models, the minimum 
value of xiim that is required for the formation of a quasi-steady Keplerian disk is given by 
xiim = Xkep = a"^/^. Let B,n,o,min denote the value of Bm,o,iim in this case. Then, Bm,o,min is 
the minimum meridional magnetic field that is required at the initial surface for Keplerian 
disk formation. When we increase the value of Xum, the pre-Keplerian region between X^^p 
and Xiim increases in size as does the radial extent of the quasi-steady Keplerian disk region 
into which it evolves. This continues until the value of xu^ reaches the value Xesc- Let 
Bm,o,opt be the value of of Bm,o,iim when xnm = Xesc = 2^/^q;~^/^. For values of xum larger 
than Xesc, although the pre-Keplerian region is larger, the radial extent of the corresponding 
quasi-steady Keplerian disk region does not increase in size because only the portion of 
the pre-Keplerian region between xn^ and Xesc evolves into the quasi-steady Keplerian disk 
region. Thus, B.rnfi,opt is smallest meridional field for which the quasi-steady Keplerian disk 
has the largest radial extent. The model with Bm,o,iim = Bm,o,opt is our optimal model. We 
denote the corresponding azimuthal field by B^^^opt- For values of Bm,o,iim > Bm,o,opt the 
radial extent of the quasi-steady Keplerian region is the same as that of the optimal model. 
Let Po,fcep(^0! ^o,fcep! 0) be the point on the initial surface that corresponds to X^ep- Then, the 
sector 'Eo,d,sk is given by 6*0,^6^ < 9o < 9o^iim. 

In applications to models, we take the initial surface, Sq, to be the sonic surface. Then, 
Vm,o = where c| = ^T/ii. Here, 3? is the gas constant and fi is the mean molecular 
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weight. Since this surface is close enough to the stellar surface we take ro ~ R. We consider 
a thin disk region with conical shock boundaries given by ^ = (7r/2) — S and 9 = (7r/2) + S 
during the fill-up stage, where S is the same for all points on the shock surface. At the end 
of the fill-up stage 6 = Sfm. Fill-up occurs on a time of order tfm, which is much longer than 
the fiow time tfiow = R/vg near the shock surface. The magnetic field emerging from the 
sector So,disA; has a dipole-type structure given by Equations (65). Since the values of 6*0,^771 
and xiim have been chosen, Equation (69) fixes the value of A. We assume that the initial 
density, po, is constant over Eq and we use the mass loss rate, M, to obtain 

Using the definition of Alfven velocity and the values of azimuthal quantities given in Sections 
4.1 and 4.2, we can write 

B^,.AU^f(M^] (73) 

and 



where f is the rotational speed at the point of the photosphere and is given by Equations 
(24) and (23) with a^, — a and 9^, — 9q. In general, Amfi and A^^ are functions of colatitude 

In order to proceed further, we need to determine either the meridional velocity, v^, 
or the density, p, of the flow as it approaches the shock boundary of the disk. When one 
of these quantities is known the other can be computed using Equation (8) in models with 
dipole-type fields. In a study where detailed properties of the wind zones are required, we 
should convert the BernouUi Equation (14) to an equation in Vm as done by Mestel (1968) or 
to an equation in p as in Sakurai (1985). Both forms of the equation arc quite complicated 
and the numerical methods that must be used to find solutions are too laborious for our 
immediate needs. Since wc require values of Vm for the fiow only near the shock surface, we 
use the /5-law for the velocity of stellar winds and assume that 

(R\^ 
1 - —J , (75) 

where /3 is a constant such that 0.5 < /3 < 2 (e.g., Lamers & Cassinelli 1999). We use 
Equation (75) to compute the flow velocity adjacent to the point X^m- Then, we compute 
the density at that point by using Equations (8) and (70). Next, we use the equations in 
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Section 10.1 to obtain the values of Am,o,iim ^-nd ^0,o,Hm, and then use Equation (71) to 
find the value of Amfi- At this stage we have the value of v4.0,o,«jm but because we have 
not specified the azimuthal field along the initial surface, we do not have any conditions on 
how Acf,fl relates to ^0,o,«jm; or how B^^q relates to -B^^o.iim- In our models, we choose the 
azimuthal field at the initial surface such that B^^q is directly proportional to v^^ at points 
on the sector So,disfc- Then, .A^^o is constant over this sector so that its value is equal to 
A(j,fl^iim- When Am,o and A^p^ are known. Equations (73) and (74) give the values of 5^,0 
and B^^o. 



13. Results 

We apply the method of solution described in Section 12 to stellar models whose basic 
properties are given in Table 1. The values of the density, and the fiow velocity at the 
initial surface are also shown in this table. We consider a variety of different cases by taking 
different combinations of values of a, ^o,/im) ^Hm, and (3. Table 2 gives the upper bounds on 
rotation rate, a, for which quasi-steady Keplerian disks can be formed when the magnetic 
field components have their optimal values. We see that the upper bound on a for Keplerian 
disk formation decreases when the terminal velocity, v^o, of the wind decreases and this can 
be explained in terms of the necessary condition (43) that we derive in Section 7, where the 
magnetic pressure in the disk plays a role even though the magnetic force is small compared 
with gravity or centrifugal force. The upper bounds on the rotation rate, a, for Keplerian 
disk formation will be larger for stars in which the normal component of the flow velocity 
into the shock surface is larger than what we have assumed in our models. 

Table 3 displays the values of the critical quantities -Bm,o,mmj Bm,o,opt and B^fl^opt foi' the 
magnetic field at the initial surface for models with selected values of a when 6*0,^771 = 65° 
and Oo^iim = 75°. We use Equation (45) to compute the fill- up times shown in the table. The 
table also shows the values of Sr-,*,L, which is the minimum field strength required by the 
meridional component of the field at the photosphere to withstand the effects of rotationally 
driven circulation. We find that the values of Bjnfi,min and Bmfl,opt are about 30 G and 70 G, 
respectively, for a BO star and about 1 G and 1.6 G for a B9 star. The magnitudes of -B<^,o,opt 
are about 15 G and 0.2 G for the same stars. Except in models of the B9 star, we find that 
the values of i?m,o,opt are less than -Br,*,L- In the B9 models, they are of the same order. In 
a rapidly rotating 03 star, -Bm,o.mm is about 500 G and -Bm,o,opt is about 1000 G. In an 06.5 
star the corresponding values are about 100 G and 200 G, respectively. In models of 03 
stars, the values of Srn,o required to form Keplerian disks are larger than 5r,*,L and in 06.5 
stars they are of the same order. At the end of the fill-up stage, the density of the disks vary 
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from about 5 x 10~^^gm cm~^ for BO stars to about 7 x 10~^^gm cm~^ for B9 stars. The 
values for 03 and 06.5 stars are, respectively, of order 10~^^gm cm~^ and 10~^^gm cm~^. 
The values of the disk density p depend on the values of M and v^o- If the mass loss rates 
near the equator for the different stellar models are higher than the values shown in Table 
1, the disk densities will be larger than the values shown in Table 3. The disk fill-up times 
have been computed for models with dfm = 10° and we find that these times are of the order 
of a few months for models of 0-type stars and about a week for B9 models. 

Table 4 shows the average values of the ratio of the magnetic force to centrifugal force 
and the ratio of the viscous force to magnetic force in the disk regions after the fill-up stage. 
Because of the increased density, the magnetic force components J?> and J-"^ in the disk 
region, given by Equations (49) and (50), become small compared with the centrifugal force 
or gravity. The results given in Table 5 for the optimal models of the B2 star show that 
when T,i,^disk is located near the latitude of 10°, its spread in latitude is only about 2° and 
when it is located near the latitude of 25° the spread is about 5°. The spread in latitude 
required for the sector Tj^^disk in optimal models is smaller when it is closer to the equator. 
Table 6 displays values of Bm,o,iim and B^fi^um for different values of Xum in models of the 
B2 star with a = 0.5 and ^o,Zim = 70°. This verifies the result we derived in in Section 12 
that the magnitudes of -8^,0 and B^^ increase when the values of are increased. 



14. Discussion 

In Section 6, we establish an important theorem connecting the angular velocity at 
points on the photosphere with the angular velocities at corresponding points in a disk 
region with no meridional motion, when the two sets of points are linked by field lines of a 
magnetic field of sufficient strength and the system is in a steady state with axial symmetry. 
This theorem is applicable in astrophysical situations in which a central rotating object 
is linked to an equatorial disk region by a magnetic field with the required strength and 
structure. In this paper, we use the theorem to discuss the formation of Keplerian disk 
regions around magnetic rotator stars. If the sector T,i,^(iisk of the photosphere rotates with 
approximately uniform angular velocity and is linked to a pre-Keplerian disk region by a 
dipole-type magnetic field of sufficient strength, then during the fill-up stage the disk region 
will have the same uniform angular velocity. In Section 7 we derive a necessary condition 
for the formation of a Keplerian disk. It requires that the normal component of the fiow 
velocity near the shock boundary must be sufficiently large for Keplerian disks to form when 
a magnetic field is present and when the rotation rate is high. This constraint results from 
the fact that magnetic pressure in the disk region plays a role in inhibiting Keplerian disk 



-31 - 



formation even when the magnetic force is small compared with gravity or centrifugal force. 

Our models require modest magnetic fields such that the meridional Alfven speeds are 
larger than the speeds of the outflow in the wind zone Wdisk between the stellar surface 
and the shock surface. Several authors (Babel & Montmerle 1997b; Donati et al 2001; 
MacGregor & Cassinelli 2003) have referred to observational evidence for the presence of 
surface magnetic fields in hot stars with strengths of several hundred gauss. We know that 
fossil fields can be present in radiative envelopes of early type stars. Spitzcr (1962) gives the 
Ohmic decay time for a stellar field with length scale L = R and appropriate resistivity to 
be 2 X 10~^^T^/^ s. If we take T = Tg//, we see that this time is much longer than the 
evolutionary time for all the stars except for the B9 star. For the B9 star they are of the 
same order. In the radiative envelopes of the stars, T is much larger than T^fj and the decay 
time will be much longer for all the stars. Recently, Charbonneau & MacGregor (2001) 
showed that magnetic fields may be generated through dynamo action at the core-envelope 
boundary of a hot star and MacGregor & Cassinelli (2003) have found that such a field 
can be transported to the surface region in a time shorter than an evolutionary time of the 
star. In view of the observational evidence and theoretical support it is reasonable for us to 
consider models with modest surface fields. 

Our approach to modeling the formation of Keplcrian disks consists of a disk fill-up 
process and a process of adjustment of super-Kcplerian disk material into Keplcrian orbits. 
Although the second process may begin during the later stages of the first process, it is 
reasonable to consider them separately because initially the density in the disk region is 
small and the magnetic field plays a role in the rotation of the disk material. Towards the 
end of the fill-up stage, the density of the disk region increases to much larger values and 
the magnetic force in the disk becomes negligible when compared with centrifugal force or 
gravity. This allows the super-Keplerian material in the pre-Keplerian disk region to move 
into Keplcrian orbits while conserving angular momentum when there is no viscous force. 
Although viscosity is negligible at the end of the fill-up stage, the onset of a weak viscous 
force during the process of diffusion of super-Keplerian material into into Keplcrian orbits 
will be helpful. 

When the meridional field S„i.o at the initial surface is only slightly larger than the 
minimum value Sm,o,mm that is required for disk formation, the Keplcrian region will be 
in the shape of a ring. For increasing values of -B„t,o the radial extent of the Keplcrian 
region increases and the region reaches its largest extent when -Bm,o equals the optimal value 
Bmfi,opt- For values of -8^,0 larger than -Brre,o,opt there is no increase in the radial extent of 
the quasi-steady Keplcrian region. When the sector Ti^^^disk of the photosphere has uniform 
rotation, the inner and outer radii of the maximal quasi-steady Keplcrian region of the disk 
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are given by a~'^/^R and 2^/^Q;'~^/^i?, respectively, where R is the stellar radius and a is 
the ratio of the rotational speed to the Keplerian speed at T,i,^cLisk- While the magnetic field 
strength at the stellar surface determines the extent of a quasi-steady Keplerian disk region, 
we see that both the inner and outer radii of a maximal disk arc given in terms of a. The 
inner radius is the distance at which there is corotation with Keplerian rotation speed. The 
outer radius of the maximal disk is determined by the point where the escape velocity is 
reached in a pre-Keplerian region and it is larger for smaller values of a. Our result for 
the maximal radial extent of a quasi-steady Keplerian region of the disk is obtained without 
including the effect of viscous forces. As shown in Table 4, the viscous force in the disk region 
is negligible at the end of the fill-up stage. If a weak viscous force plays a role during the 
transition from a pre-Keplerian to a quasi-steady Kelperian disk region, the radial extent of 
the quasi-steady region will be slightly different. There are other situations that we do not 
discuss here in which quasi-steady Keplerian disks of different radial extent can be formed. 
For example, if the magnetic field lines from the star to the disk region are confined to 
a sector S*,disA; such that the rotation speed at the innermost point of the pre-Keplerian 
disk region is super-Keplerian, then the inner radius will be larger than a~^^^R. Thus, if 
the rotation rate and the extent of a Keplerian disk region of a star can be obtained from 
observations, our models can be used to make some deductions about the magnetic field of 
the star. 

The minimum and optimal meridional magnetic fields required for the formation of 
Keplerian disks in models of B-type stars are, respectively, of order 1 G to 10 G. These fields 
are weaker than Br^^^L, which is the lower bound derived by Maheswaran & Cassinelli (1988, 
1992) for the meridional field at the stellar surface to withstand the effects of meridional 
circulation. The fill-up times in models of all the different stars that we consider in this 
paper are shorter than the corresponding circulation times. Thus, although meridional 
circulation may not affect the formation of disks in B-typc stars, its effect on the stellar field 
can infiuence some observed properties of disks. On the other hand, the situation in regard 
to 0-type stars is different. In models of 06.5 stars, meridional magnetic fields required for 
Keplerian disk formation are of the order of 100 G, which is of the same order as Br,i,,L, 
and in 03 stars they are about 1000 G, which is much larger than -B^^+^l. Although the 
fill-up times in 0-type stars arc shorter than the meridional circulation times, if these stars 
posses the meridional fields with the strengths that are required for Keplerian disk formation, 
circulation may not be involved in causing any observable changes in disks. Our numerical 
results indicate that for reasonable values of the rotation rate and surface magnetic field 
strengths, B-type stars are likely to possess Keplerian disks. This is consistent with the 
results for the distribution of Be stars obtained by Pols et al (1991). 

Our models do not require that the meridional field, Bm,o, be uniformly strong across 
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the stellar surface for the formation of a Keplerian disk. The field must have the required 

strength on the sector T^o^disk and the rotation should be approximately uniform on the sector 
^*,disk- The entire star need not be in uniform rotation. For all reasonable locations of Ti^^^isk 
in all our models, we find that its spread in latitude needs to be only a few degrees. We 
should emphasize that the assumption of axial symmetry has been made only to simplify the 
mathematics. Our model for disk formation does not depend on an aligned dipole field or 
even an axially symmetric field. Our results will qualitatively apply to oblique rotator models 
or to stars with magnetic fields consisting of fiux loops that emerge from lower latitudes and 
thread the disk region. 

The strengths of the meridional fields required for Keplerian disk formation in our 
models suggest that a possible scenario for the variation with time of properties of Keplerian 
disks in some B-type stars is the following. These stars possess magnetic fields that are 
mainly sub-photospheric and from time to time, due to magnetic buoyancy or other interior 
phenomena, meridional magnetic field lines emerge from the stellar surface. The strength 
of the meridional component of the field satisfies the condition Bm,o,min < -Bm,o < -Sr,*,L 
so that it is able to assist in the formation of Keplerian disks. After the formation of a 
disk, the meridional field lines are drawn below the stellar surface in a time scale of the 
order of the circulation time, which is longer than the disk fill-up time, tfm. This shuts 
off the supply of magnetically torqued material to the disk from non-equatorial latitudes. 
The resurgence and submersion of meridional fields across the stellar surface could explain 
some of the variations of the disk region over different time scales. We should emphasize 
that the theory that we propose for Keplerian disk formation does not depend on which 
theory of meridional circulation speeds is valid. We consider circulation only as a possible 
cause of variation of disks with time. Thus, the Be star phenomenon may be an interesting 
example where the effects of rotation, circulation and magnetic activity play a significant 
role. Obviously, stars in which i?„i.o,m.m > -Br,*.L will not exhibit this aspect of time variation. 
It is beyond the scope of the present work to consider details of the time dependent behavior 
of the surface field caused by rotation, circulation and outflow or the depletion of flux due 
to instabilities. Obviously, such a study will be greatly beneficial for a better understanding 
the role of magnetic fields in the formation, properties and variability of disks. 

Some of the basic goals of this paper are to understand how magnetic fields of rotating 
stars can influence the formation of Keplerian disks and to explain the time variability 
observed in Be stars. The MTD model of CassineUi et al (2002) has the same objectives 
and the flow regions that we consider in our models are similar to those in the MTD model. 
One of the main differences between our models and the MTD model is that we are able to 
establish a result for the angular velocity distribution during the fill-up stage of a disk. The 
MTD model uses an empirical formula for the rotational velocity in the disk region, which 
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may be appropriate in WCD models. A second difference is that the meridional fields in our 
models can be weaker than the critical strength, -Br,*,L and yet be strong enough to provide 
the torque required by the flow material to reach rotational speeds that lead to the formation 
of Keplerian disks. The MTD model requires stronger flelds because the torque is applied 
to the disk and, therefore, their strengths are larger than -Br,*,L- As we show in this paper, 
a problem with stronger fields in the disk region is that they can inhibit the formation of a 
Keplerian region. An important aspect of our models is that for the strengths of fields that 
we consider, the magnetic force in the disk region becomes negligible at the end of the fill-up 
stage. Thus, the one-armed spiral pattern of the Global Disk Oscillation model (Okazaki 
1997) can be used to explain the V/R variabihty in disks of Be stars. 

I wish to thank J. P. Cassinelli for valuable discussions and comments on the manuscript 
and gratefully acknowledge a University of Wisconsin Marathon County Foundation Summer 
Grant. 



A. Equations for Azimuthal Components and Alfven Velocities 



The value of the quantity £ for any streamline in the Northern Hemisphere is given by 
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We derive this by first substituting for C from Equation (19) in Equation (12) and then 
eliminating uj between this new equation and Equation (5). Next, we substitute for ^ and x 
from Equation (7) and use Equation (15) to replace Bm and by Am and 



We can write an expression for and in the Northern Hemisphere in the forms 
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Along an open streamline, both and B^ have removable singularities at r sin 9 — i where 
p — Pc- The values of pc and i to be inserted in these equations are given by Equations 
(18) and (Al), respectively. Expressions for Am and A^ in the Northern Hemisphere are 
obtained from Equations (8), (18), (21) and (22). These are 
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and 
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{£V(r2 sin^^)} 



l-{pc/p){iy{r^smH)} 
at points in all wind zones in the Northern Hemisphere. 
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Fig. 1. — Schematic figure of magnetic field lines, streamlines and equatorial disk during 
the fill-up stage. Material from the sector T:o,disk of the initial surface travels through the 
wind zone Wdisk and enters the disk region D. Material from So,outer flows through the wind 
zone Wouter out into space. Material from the inner sector Sumner flows into the wind zone 
W 

^ ' inner ■ 



-39- 



Table 1. Properties of the different stellar models considered 



Spectral 


M 


R 


L 




M 


Voo 


PO 


Vm.,0 


'^^,-k,eq,crit 


Type 




(Re) 




(10* °K) 


(lO-9M0yr-l) 


(km s-l) 


(gm cm~^) 


(10® cm s-i) 


(km s-l) 


03 


55 


14 


1.10 X 10® 


5.0 


9100 


3100 


1.85 X 10"" 


2.60 


592 


06.5 


29 


10 


2.31 X 10'' 


4.0 


310 


2500 


1.38 X 10-12 


2.32 


661 


BO 


15 


6.6 


4.12 X 10"' 


3.2 


27 


1300 


3.09 X 10-1=^ 


2.08 


634 


B2 


8.3 


4.5 


5110 


2.3 


0.4 


840 


1.16 X 10-1"* 


1.76 


589 


B5 


4.5 


3.5 


560 


1.5 


0.01 


580 


5.95 X 10-i'5 


1.42 


495 


B9 


2.6 


2.6 


61 


1.0 


0.0013 


460 


1.72 X 10-1® 


1.16 


437 



Note. — Entries for the basic stellar properties are from Bjorkman & Cassinelli (1993). We use L/Lq = (R/Rq)'^ {Tef f /TQ,eff)* 
to compute luminosity. 
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Table 2. Upper bounds on the rotation rate a 



Spectral 


/3 = 0.75 


(3 = 1.00 


(3 = 1.50 


Type 


a 


a 


a 


03 


0.91 


0.86 


0.72 


06.5 


0.89 


0.80 


0.65 


BO 


0.76 


0.65 


0.50 


B2 


0.65 


0.55 


0.41 


B5 


0.59 


0.49 


0.36 


B9 


0.55 


0.45 


0.33 



Note. — Values shown are for optimal mod- 
els of Keplerian disks when 9o,(im = 70°. 
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Table 3. Critical magnetic field strengths, disk fill-up times and time scales for circulation 



Spectral 


a 






B ra ,[) - rain 


Bj-,1.0, opt 


^(p,Q,opt 




^circ 


tfUl 


Pf ill ,int 


Type 




(deg) 


(deg) 


(G) 


(G) 


(G) 


(G) 


(years) 


(years) 


(grn cm~^) 


03 


0.60 


65 


70.7 


437 


893 


-63 


139 


0.53 


0.18 


9.5 X 10-11 






75 


78.5 


540 


1100 


-46 








6.2 X 10"" 


06.5 


0.60 


65 


70.7 


101 


207 


-20 


175 


0.57 


0.12 


6.0 X 10-12 






75 


78.5 


125 


256 


-15 








3.9 X 10-12 


BO 


0.60 


65 


70.7 


33 


67 


-12 


131 


0.31 


0.03 


5.0 X 10-13 






75 


78.5 


40 


83 


-9 








3.2 X 10^13 


B2 


0.50 


65 


69.6 


7.0 


13 


-2.1 


25 


1.1 


0.02 


8.7 X 10-15 






75 


77.8 


8.6 


16 


-1.6 








5.4 X 10-15 


B5 


0.45 


65 


69.1 


1.4 


2.5 


-0.40 


3.6 


6.7 


0.02 


2.7 X 10-18 






75 


77.5 


1.8 


3.1 


-0.30 








1.7 X 10-18 


B9 


0.40 


65 


68.7 


0.76 


1.3 


-0.18 


0.5 


30.6 


0.02 


7.0 X 10-" 






75 


77.3 


0.94 


1.6 


-0.13 








4.3 X lO-i'' 



Note. — We use /3 = 1.0 to compute Bmfi,min, Bm,o,opt and i?0,o,opt- Values oitfm are computed for Sfm = 10°. 
Density of the disk at the point Xj„t at the end of the fill-up stage is pint, fill 
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Table 4. Ratios of magnetic force to centrifugal force and viscous force to magnetic force 

at the end of the disk fill-up stage 



Spectral 
Type 


a 


(deg) 


ag/ ^cen 




03 


0.60 


65 


0.005 


0.056 






75 


0.005 


0.057 


06.5 


0.60 


65 


0.004 


0.045 






75 


0.004 


0.046 


BO 


0.60 


65 


0.011 


0.016 






75 


0.011 


0.015 


B2 


0.50 


65 


0.014 


0.011 






75 


0.014 


0.011 


B5 


0.45 


65 


0.017 


0.008 






75 


0.018 


0.008 


B9 


0.40 


65 


0.012 


0.010 






75 


0.013 


0.009 



Note. — We tcike v = 0.1 and /? = 1.0. Fmag / Fcen is an 
estimate of the ratio of the magnetic force to centrifugal 
force in the disk and is the average value of / {iir p f mv^) 
taken at Xf^^^ and X^^. F^iad Fmag is an estimate of 
the ratio of the viscous force to magnetic force and is the 
average value of {i-K p f mvcg / B^) at the same two points. 



Table 5. Critical magnetic field strengths, extents of optimal disk regions and latitude 

ranges of T:o,disk for the B2 star model 



do.lim = 70° a = 0.5 



Optimal model Optimal model 



a 


^kep 


^end 


^771,0, opt 










-Bfri.O,opt 












(G) 


(G) 


(G) 


(dcg) 


(dcg) 


(G) 


(G) 


(G) 


0.35 


2.01 


5.07 


22.4 


-1.4 


13.9 


60 


65.4 


12.0 


-2.4 


6.6 


0.40 


1.84 


4.64 


18.9 


-1.6 


11.3 


65 


69.6 


12.7 


-2.1 


7.0 


0.45 


1.70 


4.29 


16.2 


-1.7 


9.3 


75 


7.8 


15.8 


-1.5 


8.6 


0.50 


1.59 


4.00 


13.9 


-1.9 


7.6 


80 


81.9 


19.1 


-1.3 


10.4 


0.55 


1.49 


3.75 


12.0 


-2.0 


6.3 


85 


85.9 


26.7 


-0.9 


14.7 



Note. — 



Here /3 = 1.0 in all cases. 
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Table 6. Variation of surface magnetic field strengths and latitude ranges of T,Q^aisk with 











^kep 






(G) 


(G) 


(deg) 


1.62 


1.72 


8.1 


-1.75 


70.4 


2.02 


4.00 


13.9 


-1.85 


73.7 


3.18 


4.00 


32.8 


-2.26 


77.3 


4.77 


4.00 


61.4 


-2.86 


79.0 



Note. — Values shown are for models of the 
B2 star with a = 0.5, 6»o,am = 70° and /? = 1.0. 
In all these models, x^ep = 1-59 and Xesc = 2.00. 



